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Abstract — Several recent standards such as IEEE 802.11a/g, IEEE 802.16, and ECMA Multi- 
band Orthogonal Frequency Division Multiplexing (MB-OFDM) for high data-rate Ultra-Wideband 
(UWB), employ bit-interleaved convolutionally-coded multicarrier modulation over quasi-static fad- 
ing channels. Motivated by the lack of appropriate error rate analysis techniques for this popular type 
of system and channel model, we present two novel analytical methods for bit error rate (BER) es- 
timation of coded multicarrier systems operating over frequency-selective quasi-static channels with 
non-ideal interleaving. In the first method, the approximate performance of the system is calculated 
for each realization of the channel, which is suitable for obtaining the outage BER performance (a 
common performance measure for e.g. MB-OFDM systems). The second method assumes Rayleigh 
distributed frequency-domain subcarrier channel gains and knowledge of their correlation matrix, 
and can be used to directly obtain the average BER performance. Both methods are applicable to 
convolutionally-coded interleaved multicarrier systems employing Quadrature Amplitude Modulation 
(QAM), and are also able to account for narrowband interference (modeled as a sum of tone inter- 
ferers). To illustrate the application of the proposed analysis, both methods are used to study the 
performance of a tone-interference-impaired MB-OFDM system. 

Index terms: Error probability, channel coding, quasi-static frequency-selective fading channels, 
narrowband interference, multiband orthogonal frequency division multiplexing (OFDM), ultra wide- 
band (UWB). 
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1 Introduction 

Multicarrier communication systems based on Orthogonal Frequency Division Multiplexing (OFDM) 
have gained interest from the communications community in recent years, as evidenced by standards 
such as xDSL (digital subscriber lines), IEEE 802.11a/g (wireless local area networks) [T], IEEE 
802.16 (broadband wireless access) 0, and ECMA Multiband OFDM (MB-OFDM) for high-rate 
Ultra-Wideband (UWB) Pl E]. In general, the channel for these systems can be assumed to be 
very slowly time-varying relative to the transmission rate of the device, and can be approximated as 
quasi-static for the duration of one or more packet transmissions. As well, since most OFDM devices 
use a relatively large bandwidth, the channel is frequency selective. Simulation-based approaches to 
obtain system performance in this setting are very time consuming due to the necessity of simulating 
the system over a large number of channel realizations. 

The complexity of simulation and the wide range of practical systems operating over this type of 
channel motivate the relevance of, and interest in, analytical methods for obtaining the performance 
of coded multicarrier systems when transmitting over quasi-static frequency-selective fading chan- 
nels. There are well-known techniques for bounding the performance of convolutionally-encoded 
transmission over many types of fading channels, e.g. However, such classical bit error rate 

(BER) analysis techniques are not applicable to the OFDM systems mentioned above for several rea- 
sons. Firstly, the short-length channel-coded packet-based transmissions are non-ideally interleaved, 
which results in non-zero correlation between adjacent coded bits. Secondly, and more importantly, 
the quasi-static nature limits the number of distinct channel gains to the (relatively small) number 
of OFDM tones. This small number of distinct channel gains must not be approximated by the full 
fading distribution for a valid performance analysis, as would be the case in a fast-fading channel. 

In the quasi-static channel setting it is also often of significant interest to obtain the outage 
BER performance, i.e., the minimum expected BER performance after excluding some percentage 
of the worst-performing channel realizations |7| Section III.C-2]. Until now, one has had to resort to 
intensive numerical simulations in order to obtain the BER performance for each channel realization, 
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and hence obtain the outage BER performance |H]. 

Motivated by the considerations mentioned above, we have developed two analytical methods 
to obtain outage as well as average BER performance over ensembles of channel realizations. The 
two methods are best-suited to different types of analysis, as will be discussed in detail below. 

Furthermore, in the UWB setting, OFDM-based schemes are used as spectral underlay systems 
for previously allocated frequency bands. The resultant interference in this setting motivates the 
investigation of the effect of narrowband interference on OFDM systems. While narrowband inter- 
ference models have been considered for impulse radio and direct-sequence UWB systems there 
are no results available in the open technical literature on the effect of narrowband interference on 
coded MB-OFDM. 

Contributions: We briefly outline the main contributions of this paper: 

• We develop a method for approximating the BER of coded multicarrier systems on a per- 
realization basis ("Method I"). This method is most suitable for obtaining the outage BER, 
but can also be used to obtain the average BER performance (Section [3.3|) . 

• For quasi-static channels with correlated Rayleigh-distributed subcarrier channel gains, we 
present an alternative method ("Method M") to directly and efficiently obtain the average 
BER performance (Section \3A\f . 

• We model narrowband interference as a sum of tone interferers (a reasonable model for evalu- 
ating the effect of one or more interferers with narrow bandwidth as compared to the OFDM 
subcarrier spacing), and incorporate the effects of this interference into both analysis methods. 

• As a specific example, and motivated by the inherent interference in UWB communications, 
we provide numerical results illustrating the two BER approximation methods applied to the 
MB-OFDM UWB system. In particular, we present both outage and average BER, with and 
without tone interference (Section |H). 

• Erasure marking and decoding [TUl [TT] is a promising technique to alleviate the effects of 
interference. We study erasure marking and decoding as a mitigation technique for tone- 
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interference-impaired coded MB-OFDM (Section [4.4p . 

Related Work: In [T2], Malkamaki and Leib consider the performance of convolutional codes 
with non-ideal interleaving over block fading channels without interference. They make use of the 
generalized transfer function (GTF) [13| method in order to obtain the pairwise error probability 
(PEP). If their technique is applied to systems with a fading block length of one (equivalent to the 
quasi-static channel), their approach is similar in some ways to Method I presented in Section 
The major difference is that Method I does not require the GTF of the code, which may become 
difficult to obtain as the number of distinct channel gains (the number of blocks in the case of a 
block-fading channel) grows [13]. Instead, we apply the novel concept of error vectors, introduced 
in Section 13.11 

The PEP for uncoded and coded (across subcarrier) MB-OFDM is given in jTH]- However, the 
authors apply a non-standard UWB channel model, consider only simple codes such as repetition 
coding, and do not consider interference. 

Organization: The remainder of this paper is organized as follows. Section |2] introduces the 
OFDM transmitter and receiver models as well as the models for the channel and for the interfering 
signals. Each model is formulated quite generally, although we also mention the specific parameters 
for the MB-OFDM system, which will be the focus of numerical results presented. In Sectional we 
develop the proposed analysis methods, which allow for per-channel-realization as well as average 
error rate approximations with and without sum-of-tones interference. Analysis and simulation 
results for several practically relevant scenarios of interest for MB-OFDM are given and discussed in 
Sectional Finally, Section El concludes this paper. 

Notation: In this paper, (■)-^ and (■)^ denote transpose and Hermitian transpose, respectively, 
diag(a:;) denotes a matrix with the elements of vector x on the main diagonal, Re{-} denotes the 
real part of a complex number, E(-) denotes the expectation of a random variable, Pr{-} denotes 
the probability of some event, DFT(-) denotes the Discrete Fourier Transform, Q{-) is the Gaussian 
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Q-function and © denotes the element-wise XOR operation. Matrices and 0^ denote the 
identity matrix and the all-zero matrix of dimension r/ x 77, respectively, O^xi denotes the all-zero 
column vector of length rj, and det(-) denotes the determinant of a matrix. 

2 System Model 

In this section, we introduce the OFDM transmitter, the channel and interference models, and the 
OFDM receiver. We again note that the methods presented herein are applicable to the general 
class of coded multicarrier systems operating over quasi-static fading channels, although we have 
chosen to focus our attention on the MB-OFDM system for high data-rate UWB. 

2.1 Transmitter 

Throughout this paper we consider an A^-subcarrier OFDM system with M-ary QAM (M-QAM) 
carrying Rm = log2(M) bits per subcarrier. Figure[T]shows the relevant portions of the OFDM trans- 
mitter. The system employs a punctured convolutional code of rate Rc. We assume that the trans- 
mitter selects RcRmN random message bits for transmission, denoted by 6 = [bi 62 • • ■bR^R^N]'^- 
The vectors c and of length Lc = RmN represent the bits after encoding/puncturing and after 
interleaving, respectively. The bits are then modulated using Af-QAM on each subcarrier, and 
the resulting N modulated symbols are denoted by a? = [xi X2 ... xn]^ . 

As a specific example, we consider the MB-OFDM system proposed for the IEEE 802.15 TG3a 
high data-rate UWB standard and recently standardized by the ECMA pmHH]. The MB-OFDM 
system employs 128 subcarriers, and hops over 3 sub-bands (one hop per OFDM symbol) for first- 
generation devices. We assume without loss of generality that hopping pattern 1 of [3| is used 
(i.e. the sub-bands are hopped in order). As a result we can consider MB-OFDM as an equivalent 
384 subcarrier OFDM system. After disregarding pilot, guard, and other reserved subcarriers, we 
have N = 300 data-carrying subcarriers. 

Channel coding in the proposed standard consists of a punctured maximum free distance rate 
1/3 constraint length 7 convolutional encoder and a multi-stage block-based interleaver (see [4J for 
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details). After modulation, modulated symbols are optionally repeated in two consecutive OFDM 
symbols and/or two subcarriers within the same OFDM symbol, reducing the effective code rate 
by a factor of 2 or 4 and providing an additional spreading gain for low data rate modes. We 
can equivalently consider this repetition as a lower-rate convolutional code with repeated generator 
polynomials. The proposed standard specifies only 4-QAM with Gray labeling. We also consider 
Gray-labeled 16-QAM as a potential extension for increased data rates, and since it is also used in 
other OFDM-based standards such as IEEE 802.11a/g and IEEE 802.16. 

2.2 Channel Model 

We will assume that the OFDM system is designed such that the cyclic prefix is longer than the 
channel impulse response. Thus, we can equivalently consider the channel in the frequency domain, 
and denote the subcarrier gains by h = [hi h2 ... hj^i]. We also include the frequency-domain 
interference signal J (see Section The transmitted symbols x pass through the fading channel 
H = diag(/i), and the length- vector of received symbols r (after the DFT) is given by 

r = Hx + J + n , (1) 

where n is a vector of independent complex additive white Gaussian noise (AWGN) variables with 
variance Aq. We denote the energy per modulated symbol by Eg = RcRmEb, where Eh is the 
energy per information bit. 

For a meaningful performance analysis of the MB-OFDM proposal, we consider the channel 
impulse response model developed under IEEE 802.15 for UWB systems [1^1 (a version of the Saleh- 
Valenzuela model [T7] modified to fit the properties of measured UWB channels). Multipath rays 
arrive in clusters, with exponentially distributed cluster and ray interarrival times. Both clusters 
and rays have decay factors chosen to meet a given power decay profile. The ray amplitudes are 
modeled as lognormal random variables, and each cluster of rays also undergoes a lognormal fading. 
To provide a fair system comparison, the total multipath energy is normalized to unity. As well, the 
entire impulse response undergoes an "outer" lognormal shadowing. The channel impulse response 
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is assumed time invariant during the transmission period of (at least) one packet. We consider the 
UWB channel parameter set CMl for short-range line-of-sight channels [T^. We also note that, due 
to the relative complexity of this channel model, even BER analysis for uncoded UWB systems has 
proven challenging [TH] , 

The outer lognormal shadowing term mentioned above is irrelevant for the fading characteristics 
as it affects all tones equally. Denoting the lognormal term by G, we obtain the corresponding 
normalized frequency-domain fading coefficients as 

/i" = h/G . (2) 

The elements of /i" are well-approximated as zero-mean complex Gaussian random variables ^\ 
EH]. This allows us to apply analysis assuming correlated Rayleigh fading coefficients (see Section lT4|l 
to the UWB channel without lognormal shadowing, and then average over the lognormal shadowing 
distribution in order to obtain the final system performance over the UWB channel. We note that this 
is only relevant for the method in Section [TH — for the realization-based method (see Section IT5|) 
the distribution of h is not important. 

2.3 Interference Model 

We model narrowband interference as the sum of A^^ tone interferers 

^{t) = Y,^k{t) , (3) 

k=l 

where the equivalent complex baseband representation of the A;th tone interferer with amplitude a^, 
frequency fk, and initial phase (pk is given by 

tk{t) = ake^^'-f^'+^-^ . (4) 

Assuming that the interference i{t) falls completely in the passband of the receiver filter before 
sampling, we form the discrete-time equivalent interference by sampling i(t) with the OFDM system 
sampling period T, and obtain (for one OFDM symbol) the sample vector 

i=m z{T) z{2T) ... z{{N-l)T)f. (5) 
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Therefore, the frequency-domain equivalent J of the interfering signal considered in ((T]) is given 

by 

J = DFT(*) . (6) 

We note that, due to the finite-length DFT window, each single-tone interferer is convolved by 
a sinc-function in the frequency domain. If fk is equal to one of the subcarrier frequencies, only one 
subcarrier is impaired by the interferer ik{t) (since the interferer will be zero at the other subcarrier 
frequencies). On the other hand, if fk happens to lie between two subcarriers, the tone interferer 
will affect several adjacent subcarriers. 

2.4 Receiver 

The relevant portions of the OFDM receiver are shown in Figure [T] We assume perfect timing and 
frequency synchronization. The receiver employs a soft-output detector followed by a deinterleaver 
and a depuncturer. After possible erasure marking based on knowledge of /jt, 1 < /c < Ni (see 
Section for details), standard Viterbi decoding results in an estimate b = [bi 62 ... hR^RmNY 
of the original transmitted information bits. 

3 Performance Analysis 

in this section, we present two methods for approximating the performance of coded multicarrier 
systems operating over frequency-selective, quasi-static fading channels and impaired by sum-of- 
tones interference. The first method (Section 13. 3p is based on approximating the performance of 
the system for individual channel realizations. The main strength of this method is that it can be 
used to obtain the outage BER performance (the standard performance measure considered in MB- 
OFDM systems PJEIIH]). While the first method can also be used to obtain the average BER over 
an ensemble of channel realizations, the second method (Section f3.4|) . which is based on knowledge 
of the correlation matrix of the frequency-domain channel gains, can be used to directly obtain the 
average performance without the need to consider a large ensemble of channels. Both methods are 
based on considering the set of error vectors, introduced in Section 13.11 and the PEP of an error 



Snow et al.: Error Rate Analysis for Coded Multicarrier Systems over Quasi-Static Fading Channels 8 



vector, given in Section 13.21 

One major problem in the analysis of M-QAM schemes with M > 4 is that the probability of 
error for a given bit depends on the whole transmitted symbol (i.e., it also depends on the other bits 
in the symbol). For this reason, for the combination of convolutional coding and M-QAM it is not 
sufficient to adopt the classical approach of considering deviations from the all-zero codeword only. 
In theory, one must average over all possible choices for c. Since this is computationally intractable, 
we simply assume the transmitted information bits b (and hence x) are chosen randomly. For 
M = 4 (where the joint linearity of code and modulator is maintained) this is exactly equivalent 
to considering an all-zero codeword. In the case of M > 4, we have verified for various example 
scenarios that, for the two analysis methods proposed below, a random choice of b well-approximates 
the true system performance. 

3.1 Error Vectors 

Let £ be the set of all L vectors ei {1 < i < L) of code output (after puncturing) associated with 
input sequences whose Hamming weight is less than u^max. i e-. £ = {^1,^2, ■■■ ^^l}- We limit 
£ to only simple vectors ^Ij, where the error path contains only one deviation from (and return 
to) the all-zero trellis state. Let le and ae be the length of and the number of information bit 
errors associated with e^, respectively. As with standard union-bound techniques for convolutional 
codes the low-weight terms will dominate the error probability. Hence, it is sufficient to choose a 
small Wmax — for example, the punctured MB-OFDM code of rate Rc = 1/2 [4j has a free distance 
of 9, and choosing uJmax = 14 (resulting in a set of L = 242 error vectors of maximum length 
/ = 60) provides results virtually identical to those obtained using larger Wmax values. 

We term an "error vector" and £ the set of error vectors. Note that can be straightforwardly 
obtained from the transfer function of the code, without resorting to the GTF [13J approach adopted 
in jT2], and is also independent of the number of distinct channel gains (or the number of blocks in 
the context of [12j). 
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3.2 Pairwise Error Probability (PEP) for an Error Vector 

We consider error events starting in a given position i of the codeword {I < i < Lc). For a specific 
error vector (1 < £ < L), form the full error codeword 

q.^^ = [0 ^ ■ ■ ^ p ^ . ■ OJ ^ (7) 

i—l li Lc—lg—i+l 

of length Lc by padding with zeros on both sides as indicated above. Given the error codeword 
qfj^ and given that codeword c is transmitted, the competing codeword is given by 

Vi^l = C © Qi^^ . (8) 

The decoder employs a standard Euclidean distance metric (i.e., the interference signal is assumed 
to be unknown for calculation of the metric). Letting be the vector of M-QAM symbols 
associated with vj^. (the interleaved version of Vn), and recalling that x is the modulated symbol 
vector corresponding to the original codeword c, the PEP for the ith error vector starting in the ith 
position, i.e., the probability that Vi^e is detected given that c was transmitted, is given by 

PEPi,£(/f , J) = FT{\\r - HxW^ > \ \r - Hzi,i\\^ \H,J} . (9) 

In Sections 13.31 and we will obtain various forms for this general expression. 

3.3 Per-realization Performance Analysis ("Method 1") 

In this section, we obtain an approximation of the BER for a particular channel realization H = 
diag(h,) and interference J, which we denote as P{H, J). For simplicity, we refer to this method as 
"Method I" in the remainder of this paper. As noted above and discussed in more detail in Section 
13.3.31 the main strength of this method is the ability to obtain the outage BER of coded OFDM 
systems. 

3.3.1 Pairwise Error Probability (PEP) 

The PEP for an error vector (1 < £ < L) with the error event starting in a position i {1 < i < Lc) 

is given by (j^. For a given H and J, and after some straightforward manipulations, we obtain the 
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PEP,,,(/f , J) = g Z' ^" ^ = ^ . (10) 



expression 



It is insightful to examine ()10p for two special cases: 

• Ao (the low-noise region): In this case, there are two possible outcomes. If the numerator 
in (II Oil is positive, we have the Q-function of a large positive value and thus PEP 0. 
However, if the interference J causes the numerator to become negative, we have the Q- 
function of a large negative value and thus PEP — > 1. That is, we either (depending on J) 

will surely make an error, or will surely not make an error. 

• J = Onxi (no interference): Here we can simplify (II Op to obtain 




PEP,,(if, J) = g (^y '^^%Ar^"'^^' I • (^^^ 

3.3.2 Per-realization BER 

The corresponding bit error rate for the ith error vector, starting in the ith position, is given by 

Pi^e{H,J) = a,-FEF,^,{H,J) . (12) 

Summing over all L error vectors, we obtain an approximation of the BER for the ith starting 
position as 

L 

P,{H,J) = Y,P^AH.J) . (13) 

£=1 

We note that ()13|) can be seen as a standard truncated union bound for convolutional codes (i.e., it 
is a sum over all error events of Hamming weight less than cUmax)- We also note that we can tighten 
this bound by limiting Pi to a maximum value of 1/2 before averaging over starting positions 
Finally, since all starting positions are equally likely, the BER P{H,J) can be written as 



i=i L e=i 

Table [U contains pseudocode to calculate P(H,J) according to (|T4|l 



(14) 
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3.3.3 Average and Outage BER 

The average BER for a given interference can be obtained by averaging ()14j) over a (large) number 

Nc of channel realizations, where the ith channel realization is denoted by {1 < i < N^), as 



As mentioned previously, Method I also readily lends itself to the consideration of the outage 
BER, a common measure of performance for packet-based systems operating in quasi-static channels 
|7]. The outage BER provides a measure of the minimum performance that can be expected of the 
system given a specified XVo outage rate, and is often employed in UWB system performance studies 
We evaluate ()14p for a set of channel realizations H = {Hi,l < i < Nc}. The worst- 
performing XVo of realizations are considered in outage, and those channel realizations are denoted 
by TYout- Denoting the remaining (100 — X)% of channel realizations by Hin, the outage BER is 
given by 



In Sectional we will focus on results for fixed values of signal-to-interference ratio (SIR), interferer 
amplitude a^, and interferer frequency f^. However, in order to remove the effect of the interferer 
initial phase, we will average (II 5p and ()16p over 32 uniformly-distributed values of 0^ G [0,27r). 

3.4 Average Performance Analysis ("Method M") 

In this section, we propose a method, based on knowledge of the frequency-domain channel corre- 
lation matrix, which can be used directly in order to obtain the average BER performance of coded 
multicarrier systems. The advantage of this method is that it allows for simple and direct evaluation 
of the average BER, without the need to evaluate the BER of many different channel realizations 
as in Method I, cf. ()15j) . For simplicity, we refer to this method as "Method M" in the remainder of 
this paper. 

For this method we will explicitly assume that the elements of h are Rayleigh-distributed and have 
known correlation matrix (in practice, S/i^j can be obtained from actual channel measurements, 




(15) 



max P{Hi, J) . 



(16) 
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or can be numerically estimated by measuring the correlation over many realizations of a given 
channel model). As noted in Section IT^ the channel models for OFDM-based UWB communication 
satisfy this assumption. 
3.4.1 Average PEP 

Noting that only the r]i/ non-zero terms of {x — Zi^i) in (j^ contribute to the PEP (and suppressing 
the dependence of 77 on i and i for notational clarity), we let x' , z'-^, H' = diag(h'), J', and n' 
represent the transmitted symbols, received symbols, channel gains, interferences, and AWGN noises 
corresponding to the 77 non-zero entries of (a; — Zi^i), respectively, and form Hh'h' by extracting the 
the elements from S^.^ which correspond to h' . Letting D = diag{x' — z[f) be the diagonal matrix 
of non-zero entries and g = H\x' — z'-^) = Dh! , we have 







^99 ~ DTih'h'D^ 



(17) 
(18) 



i.e., the distribution of g is zero-mean complex Gaussian with covariance matrix R, 



99- 



We would like to obtain the average PEPj^£ for the £th error vector, starting in the ith position. 
Rewriting including only the contributing terms, we obtain 



PEP,, 



Ft {\\r' - H'z'Jl^ - \\r' - H'x'W" < 0} , 
Pr {gg'' - g{J' + n' f - (J' + n')g^ < 0} 
Pr{A,,,(D) <0} , 



(19) 



(20) 



where Ai^i{D) = Ay and 



y 



g 

J' + n' 



We adopt the Laplace transform approach 
two typical narrowband channel situations: 



A 



'I Tj Ojy 



to determine Pr {Aj ^(Z)) < 0}, and consider 
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Case 1 — ak constant (non-faded interferers): In this case, we note y has mean fiyy and 
covariance matrix Ryy, which are given by 



Hyy = E{y) 



J' 



R 



yy 



Rgg Or, 

0^ 7Vo7^ 



The Laplace transform of Ai^£{D) is given by 

eXp[-S^l^y{A~^ + SRyy^^tlyy] 



(21) 



det(/2r, + sRyyA) 

Case 2 — at independent Rayleigh faded interferers: In this case, ]E(y) = 02,,xi, and we have 



R 



yy 



Rgg 0^ 



0,, Rj'J' + -Afolri 



(22) 



where Rj'j' =¥. (J J ), and the Laplace transform of Aj£(Z>) is given by 

1 



det(/2r, + sRyyA) 

In either case, the average PEP for the fth error vector starting in the ith position is given by 



(23) 



PEP,,, = Pr {A,,,(D) < 0} 



C+JOO 

^iAs) — 



(24) 



C-JOO 



where c is in the convergence region of (^u^s). We note that may be solved efficiently either 
analytically or via numerical integration using a Gauss-Chebyshev quadrature rule 

3.4.2 Average PEP without Interference 

An alternative form for the average PEP can be obtained for the special case of J = Oatxi- From 

(|TT|) . and following Eq. (7)], we can write the average PEP for the £th error vector starting in 
the zth position as 



PEP 



1 



7r/2 



det ( -M^^ + L 



d9 



(25) 





It can be shown that the Laplace transform approach with J = Onxi leads to an equivalent 
expression. 
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3.4.3 Average BER 

Given the average PEP according to either l^^ or ()25p . the corresponding bit error rate for the ith 
error vector, starting in the ith position, is given by 



Pi,e = ae ■ PEP,,, . (26) 
Summing over all L error vectors, the BER for the ith starting position can be written as 

L 

P^ = Y,h,,. (27) 

£=1 

Finally, since all starting positions are equally likely to be used, the average BER P can be written 
as 

^-^Efl-iEEfl.. (28) 

^ i=l ^ 1=1 1=1 

Table[T]contains pseudocode to calculate P according to ()28|) . Note that, since Pj ^ in ()28p is already 
averaged over H, we cannot upper-bound it by 1/2 as we did in (I14p for Method I. This implies 
that the result for Method II may be somewhat looser than that for Method I (see also Section 

4 Numerical Results 

in this section, we present numerical results for Methods I and II introduced above. We focus on the 
particular case of MB-OFDM operating in the CMl UWB channel. As mentioned in Section 12.21 
for Method II we include the effect of "outer" lognormal shadowing by numerically integrating the 
results of ()28p over the appropriate lognormal distribution | jl6| . 

4.1 No Interference 

In Figure |21 we present the 10% outage BER as a function of Eij/Mq obtained using Method I 
(lines), as well as simulation results (markers) for different code rates and modulation schemes using 
a set of 100 UWB CMl channel realizations with lognormal shadowing, where Eb denotes the mean 
received energy per information bit over the ensemble of channels. The 10% outage BER is a 
common performance measure in UWB systems, cf. e.g. t4|. We can see that Method I is able to 
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accurately predict the outage BER for 4-QAM and 16-QAM modulation schemes and a variety of 
different code rates, with a maximum error of less than 0.5 dB. It is also important to note that 
obtaining the Method I result requires significantly less computation than is required to obtain the 
simulation results for all 100 UWB channel realizations. For example, it took about 15 minutes to 
obtain one of the analytical curves of Figure |21 while it took approximately 48 hours to obtain the 
corresponding simulation results on the same computer. 

Figure El illustrates the average BER as a function of Eb/Mo foi' 4-QAM and 16-QAM with code 
rates Rc = 1/2 and 3/4 using two approaches: Method I with an average over 10,000 channel 
realizations (dashed lines), and the direct average from Method II (solid lines). As expected, the 
two methods are in close agreement at low BER. The deviation between the two results at higher 
BER is due to (a) the loosening effect of the averaging of Method II over the lognormal distribution, 
and (b) the fact that Method I is somewhat tighter due to the upper-bounding by 1/2 in (I14|) . 

A Caution to System Designers: We should note that 100 channel realizations (standard for 
MB-OFDM performance analysis j8|) may not be sufficient to accurately capture the true system 
performance. Figure El (solid lines) shows the average BER with respect to E^/Mq for four different 
sets of 100 UWB CMl channel realizations, obtained via Method I. For comparison, the average 
performance obtained via Method II is also shown (bold solid line). We can see that the average 
system performance obtained using sets of only 100 channel realizations depends greatly on the 
specific realizations which are chosen. Similarly, Figure El illustrates the 10% outage BER with 
respect to Eb/J\fo for four different sets of 100 UWB CMl channel realizations, obtained via Method 
I (dashed lines). For comparison the 10% outage BER obtained using a set of 1000 realizations is 
also shown (bold dashed line). We see that the outage BER curves, while less variable than the 
average BER curves, are still quite dependent on the selected channel realization set. 

Based on the results above, it seems that performance evaluation for systems operating in quasi- 
static channels using only small numbers of channel realizations may be prone to inaccurate results. 
This is one of the main strengths of the two methods presented in Section El the performance can 
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easily be evaluated over any number of channel realizations (Method I), or the average performance 
can be directly obtained (Method II), without resorting to lengthy simulations. 

4.2 Non- Faded Tone Interference 

In this section, we present results for non-faded tone interference, specifically focusing on the MB- 
OFDM system operating at 320 Mbps (i?c =1/2 after puncturing) with 4-QAM modulation over 
the CMl channel. We concentrate on the case of A'^j = 1 interferer, in order to examine the effect 
of the interferer frequency fi and the signal-to-interference ratio^ 



Without loss of generality we place /i between the 52°"^ and 53'''^ MB-OFDM subcarriers. 

In Figure El we consider E^/Mq = 17 dB, SIR = 19 dB, and focus on the effect of varying 
/i. We show the average BER for five different sets of 100 channel realizations (dashed lines), 
obtained via Method I. The markers (*) indicate simulation results which correspond to, and are in 
good agreement with, the Method I results for the set of 100 channels indicated by a bold dashed 
line. Figure El indicates that the best-case performance is obtained when fi lies exactly between 
two OFDM subcarriers (interferer position 52.5), and the performance degrades as /i approaches 
a subcarrier frequency. We also note that, as seen in Figure El for the no-interference case, the 
performance obtained using Method I can be quite variable when considering small sets of channel 
realizations. 

Figure El shows the BER versus -E^/A/q for one non-faded interferer at positions 52.0 (solid 
lines) and 52.5 (dashed lines) with SIR = {28, 23, 21, 19, 17, 15} dB, obtained using Method II. For 

comparison, the no-interference (SIR = oo) performance from Method II (bold solid line) is also 

^Note that the SIR according to this definition is an average over all the subcarriers, so the SIR for a specific 
subcarrier may be much higher/lower than the average. For example, in the 384-subcarrier MB-OFDM system with 
one interferer directly on a subcarrier, the SIR of the affected subcarrier will be « 26 dB lower than the average SIR 
(since the interference on all other subcarriers is zero). 



SIR 




(29) 
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shown. This figure clearly illustrates the performance degradation associated with decreasing SIR. 
As seen in Figure El the best-case performance is obtained when fi lies exactly between two OFDM 
subcarriers, while the performance degrades as /i approaches a subcarrier frequency. 

4.3 Rayleigh-faded Tone Interference 

We now consider the effect of A'j = 1 Rayleigh-faded interferer, in order to compare the relative 
effects of interference with those of the non-faded interferer in Section 14.21 Figure |71 shows the 
BER versus E^/N'o obtained using Method II for the same interferer positions and SIR values as 
in Figure El By comparing Figures Eland |71 we can clearly see that Rayleigh-faded tone interferers 
have a larger effect on the BER performance than non-faded tone interferers. For example, at 
Eb/J^o = 17 dB, SIR = 21 dB and interferer position 52.5, the BER with one non-faded interferer 
is approximately 10"^, while the BER with one Rayleigh interferer is approximately 2.3x10""^. Even 
for relatively high SIR = 28 dB, one Rayleigh tone interferer at position 52.0 causes a much larger 
effect than the non-faded tone interferer at the same SIR. 

4.4 Interference Mitigation by Erasure Marking and Decoding 

In OFDM systems where interference impacts a small number of subcarriers, one simple and practical 
method of interference mitigation is to erase the information bits carried on the most-affected 
subcarriers (proposed in e.g. pTO], as well as more advanced joint marking and decoding in [ill]). 
In order to study the potential performance of such an erasure technique, we consider the use of 
a genie which erases the subcarriers with largest interference powers. In the framework of analysis 
of Sectional subcarrier erasures can be considered as additional puncturing and easily incorporated 
into both analysis methods. 

Figure IHl illustrates the average BER versus Eb/Mo for 0, 1, and 2 subcarrier erasures, obtained 
using Method II. One non-faded interferer is placed at positions 52.25 (dashed lines) and 52.5 
(solid lines), with SIR of 15 and 19 dB. As can be seen from this figure, the use of a small 
number of subcarrier erasures rapidly decreases the effect of the tone interference and allows the 
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interference-impaired system performance to approach the no-interference performance (bold solid 
line). Focusing on the case of position 52.5 (solid lines), we can see that using only one erasure has 
a small effect on the resultant BER. This is due to the windowing effect of the DFT at the OFDM 
receiver (see Section fl^ . which results in interfering signal power being symmetrically distributed 
amongst a number of subcarriers. However, once the two largest equal-interference-power subcarriers 
are erased, performance improves dramatically. On the other hand, when the tone interferer is at 
position 52.25, a large portion of the interference power is in one subcarrier, so even one erasure can 
provide a substantial performance improvement. We should also note that if the interferer happens 
to be exactly at the subcarrier frequency, one subcarrier erasure will suffice to totally remove the 
effect of the interference. 

We conclude by returning once again to the consideration of outage BER obtained via Method I. 
In FigureEl we consider one non-faded interferer at position 52.5, and show the number of subcarrier 
erasures required to maintain the 10% outage BER < 10^^ for varying SIR and different values of 
Eh/Mo- As expected, decreasing the SIR results in a higher required number of erasures to maintain 
the target BER. Unfortunately, a large number of erasures compromise the code's error correcting 
capability. As can be seen from Figure^l eventually too many erasures weaken the code sufficiently 
that, even with the effects of interference mostly removed, the code is not able to maintain the 
required target BER. Figure El also shows that providing an increased SNR margin allows the MB- 
OFDM system to compensate for a larger amount of interference. 

5 Conclusions 

In this paper, we have presented two methods for evaluating the performance of convolutionally- 
coded multicarrier systems employing M-QAM, operating over frequency-selective, quasi-static, 
non-ideally interleaved fading channels and impaired by sum-of-tones interference. The realization- 
based method ("Method I") presented in Section IT^ estimates the system performance for each 
realization of the channel with an arbitrary fading distribution, and is suitable for evaluating the 
outage performance of systems. The method presented in Section ITM ("Method M"), based on 
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knowledge of the correlation matrix of the Rayleigh-distributed frequency-domain channel gains, 
allows for direct calculation of the average system performance over the ensemble of quasi-static 
fading channel realizations. 

These two novel methods allow for analytical evaluation of the performance of a class of systems 
of considerable practical interest, whose performance evaluation was previously only possible via 
intensive numerical simulations. The results in Section El demonstrate that the proposed methods of 
analysis provide an accurate measure of the system performance and allow for much greater flexibility 
than simulation-based approaches. 

We have also shown that the MB-OFDM system (and OFDM systems in general) may be 
significantly impacted by the effect of tone interference, but that this performance degradation can 
be mitigated to a large extent by the use of erasure marking and decoding at the receiver, provided 
that the receiver can obtain knowledge of which subcarriers are impaired by the interferers. 
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Tables and Figures: 



Table 1: Pseudocode for the two analysis methods. 



Method 1 


Final BER is P (for given H, J). 


1 


P := 


2 


for i := 1 to Lc do 


3 


Pi := 


4 


for i := 1 to L do 


5 


form q^^ as per (0) 


6 


form as per <^ 


7 


form vj^ and Zi^£ from Vi^£ 


8 


calculate PEPj^ as per (flOl) 


9 


calculate Pij as per (11 2D 


10 


-Pi := -Pj + -Pi,^ 


11 


endf or 


12 


P ■= P + min(i,Pi) 


13 


endf or 


14 


P := P / 


Method II 


Final BER is P. 


1 


P := 


2 


for i := 1 to Lc do 


3 


for £ := 1 to L do 


4 


form q^^ as per ([Tl) 


5 


form Vi^i as per (IHl) 


6 


form vj^ and from 


7 


form x'i,^, h,', J' and S^./^,/ 


8 


compute Z) := diag(a3' — z'i^i) 


9 


compute g = Dh' and Rgg := DHh/h' 


10 


form and/or Rj'j' (as required) 


11 


form Ryy and A 


12 


form $i,£(s) as per either (l2TD or (l23l) 


13 


calculate PEPi,^ as per (1^ 


14 


calculate Pj^ as per (I26D 


15 


P := P + Pi,i 


16 


endf or 


17 


endf or 


18 


P := P / 



Snow et al.: Error Rate Analysis for Coded Multicarrier Systems over Quasi-Static Fading Channels 22 



b 


Encoder / 


c 


Interleaver 




A/-QAM 


X 




Puncturer 






Mapper 





Transmitter 



r 


A/-QAM 
Demapper 




Deiiiteiieaver 




Depuncturer / 
Decoder 


6 













Receiver 



Figure 1: Relevant portions of the OFDM transmission system. 
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Figure 2: 10% outage BER vs. Eb/Mo ^fom Method I (lines) and simulation results (markers) for 
different code rates and modulation schemes. UWB CMl channel. Code rates 1/4 and 1/8 include 
repetition. No interference (J = Oivxi)- 
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Figure 3: Average BER versus Et/Mo for 4-QAM and 16-QAM with code rates Rc = l/2 and 3/4. 
Solid lines: Direct average from Method II. Dashed lines: Method I with an average over 10,000 
channel realizations. UWB CMl channel. No interference (J = Oatxi)- 



Snow et al.: Error Rate Analysis for Coded Multicarrier Systems over Quasi-Static Fading Channels 25 




Figure 4: Average BER (solid lines) and 10% outage BER (dashed lines) versus Eh/Mo for four 
different sets of 100 channels using Method I. For comparison: average BER from Method II (bold 
solid line), and 10% outage BER for a set of 1000 channels (bold dashed line). UWB CMl channel, 
Rc = 1/2, 16-QAM. No interference (J = O^vxi)- 
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Figure 5: Average BER versus interferer position for £'f,/A/o = 17 dB and SIR = 19 dB with 
one non-faded interferer. Results shown: Method I with different sets of 100 channel realizations 
averaged over 32 phases 0i G [0, 27r) (dashed lines); simulation results for 100 channel realizations 
corresponding to bold dashed line (markers). UWB CMl channel, = 1/2, 4-QAM. 
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Figure 6: Average BER versus E^/Mq for various SIR with one non-faded interferer, obtained using 
Method II. Interferer positions 52.0 (solid lines) and 52.5 (dashed lines). For comparison: SIR = oo 
(J = Otvxi) from Method II (bold solid line). UWB CMl channel, Rc = 1/2, 4-QAM. 
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Figure 7: Average BER versus E^/Mq for various SIR with one Rayleigh-faded interferer, obtained 
using Method II. Interferer positions 52.0 (solid lines) and 52.5 (dashed lines). For comparison: 
SIR = oo (J = Ojvxi) from Method II (bold solid line). UWB CMl channel, = 1/2, 4-QAM. 
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Figure 8: Average BER versus Eb/Mo for {0,1,2} subcarrier erasures. One non-faded interferer, 
positions 52.25 (dashed lines) and 52.5 (dashed lines) and SIR = {15,19} dB, obtained using 
Method II. For comparison: SIR = oo (J = Oatxi) from Method II (bold solid line). UWB CMl 
channel, Rc = 1/2, 4-QAM. 
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Figure 9: Required number of subcarrier erasures to maintain 10% outage BER < 10~^ for various 
SIR and E^/Mq. One non-faded interferer, position 52.5, average over 32 phases (pi E [0, 27r), 
obtained using Method I with 1000 UWB CMl channel realizations, Rc = 1/2, 4-QAM. 



